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ABSTRACT 

We present new constraints on the merger history of the most X-ray luminous cluster of 
galaxies, RXJ1347. 5-1145, based its unique multiwavelength morphology. Our X-ray analysis 
confirms the core gas is undergoing “sloshing” resulting from a prior, large scale, gravitational 
perturbation. In combination with extensive multiwavelength observations, the sloshing gas 
points to the primary and secondary clusters having had at least two prior strong gravitational 
interactions. The evidence supports a model in which the secondary subcluster with mass 
M=4.8±2.4x 10 14 Mq has previously (>0.6 Gyr ago) passed by the primary cluster, and has 
now returned for a subsequent crossing where the subcluster’s gas has been completely stripped 
from its dark matter halo. RXJ1347 is a prime example of how core gas sloshing may be used 
to constrain the merger histories of galaxy clusters through multiwavelength analyses. 
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1. Introduction 


Clusters of galaxies represent the end state 
in the hierarchical growth of the largest gravita- 
tionally bound structures in the Universe. Out- 
side their cores, the (near) universal temperature 
and mass profiles of relaxed clusters make them 
ideal tools for constraining cosmological parame- 
ters (e.g. wq and erg) whic h govern the growth o f 
structure in the universe ( Vikhlinin et al • 2009). 
We see this growth of structure in la rge sc ale arn- 


ulations (e.g. Millennium I and I I; ISpringel et al 


2005 : Bovlan-Kolchin et al. 20091 ). however, there 
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are relatively few corresponding observational con- 
straints in place. This is particularly true for 
matching the predicted rates at which massive 
clusters assemble with direct observational con- 
straints on those rates. For individual clusters, 
one should use empirically derived initial condi- 
tions in order to constrain their merger histo- 
ries and attempt to reproduce the observations 
in h and, as was done rece ntly for the Virgo clus- 
ter ( Roedieer et al.1 201l[) . One such cluster for 
which there exists a wealth of observational data 
is the most luminous X-ray cluster, RXJ1347.5- 
1145 (z=0.451). 

High resolution hydrodynamic simulations pre- 
dict the effects that multiple gravitational inter- 
actions will have on galaxy clusters. In particu- 
lar, the resulting morp holo gy o f t he l o w en trop-y 


lar, tne resulting morpnol ogy o t tiie l ow en tropy 
cluster core gas J Ascasibar_&Maxkevitc^ 2006^, 


Po ole et aDl2008t IZuHone et all 2 01^ Roedieer et al 


2 0111) . In this paper, we use the extensive multi- 
wavelength observations, along with simulations, 
to propose a merger history for RXJ1347 which 
predicts the current morphology of the cluster in 
optical, X-ray and millimeter Sunyaev-Zeldovich 
(SZ) observations. In the following section, we 
describe the phenomenon of core gas sloshing in 
galaxy cluster cores. In fj2j we present a brief dis- 
cussion of the extensive multiwavelength coverage 
and analysis of the cluster. In )}3] we discuss our 
analysis of the Chandra X-ray observations and 
the identification of the sloshing gas in the cluster 
core. Finally, in jj4] we present new constraints 
on the merger history of the cluster based on the 
accumulated multiwavelength observations. 

1.1. Core Gas Sloshing 

The model for core gas sloshing has been ex- 

plored numerically in several work s (Tittlev & Henriksen 


2005:lAscasibar fc Markev 


tchl2006 : ZuHone et al 


ject (the subcluster) falls in on either an elliptical 
or hyperbolic orbit (the exact orbit is irrelevant 
for the subcluster’s first crossing). As the subclus- 
ter falls in, the primary cluster (its dark matter 
halo and cluster gas) experience a gravitational 
force towards the incoming object. The subclus- 
ter reaches pericenter in its orbit, at which point 
the gravitational force is at its peak, and the pri- 
mary cluster reaches its maximum displacement 
from its initial position. 

The wake of the passing subcluster directs the 
ambient velocity field towards the primary clus- 
ter, while at the same time the primary clus- 
ter continues to move towards the subcluster. 
This leads to ram pressure on the primary clus- 
ter gas, causing it to compress along the lead- 
ing edge of its trajectory. As the subcluster 
recedes, the velocity field around the primary 
cluster changes rapidly, removing the ram pres- 
sur e, a nd resultin g in a ram pressure “slingshot” 

( Hallman fe Markevitchll2004 ) , where the primary 


cluster gas moves quickly ahead of the potential 
minimum, causing a separation between the two. 
It is this initial separation that leads to the sub- 
sequent production of large scale asymmetries in 
the gas distribution. 

After its initial separation, the primary clus- 
ter gas density peak (its “cool core”) falls back 
towards the potential minimum, which has itself 
been displaced from its initial position. The den- 
sity peak moves against the flow of its own trail- 
ing gas, resulting in ram pressure along its lead- 
ing edge. This ram pressure causes a significant 
(> 2x) compression of the gas along the edge, 
revealing itself as a crescent-shaped edge in the 
X-ray surface brightness distribution of the core 
gas. Upon reaching the apocenter of its oscilla- 
tion about the potential minimum, the cool core 
reverses direction and again falls back towards 


vationally in many galaxy groups (Randall et al. own trailing gas, resulting in the formation of an- 

2009; 

Machacek et al. 

20101) and clusters (IMarkevitch et al.^he 1 ' ram pressure front. These fronts have been 

2001; 

Churazov et al. 

2003: Markevitch Hr. Vikhlinin dubbed “cold fronts” ( V ikhlinin et al. 2001), as 

2007; 

Owers et al. 2009| Johnson et al. 201(1 and opposed to shock fronts, because they delineate 


Johnson et al. in preparation). 

The sloshing phenomenon begins when two 
group or cluster sized halos cross nearby one an- 
other (<1 Mpc depending on their masses). As 
seen from the reference frame of the more massive 
object (the primary cluster), the less massive ob- 


the edge of the cool core’s oscillatory motion, with 
the coolest gas interior to the front and the warmer 
gas outside. This oscillation and front production 
continues for long periods (>1 Gyr) in the core 
but is progressively damped with each crossing of 
the density peak through the potential minimum. 


2 


The above describes the essential components 
of the radial motion of the density peak about 
the potential minimum. However, because angu- 
lar momentum is transferred from the subcluster 
in this interaction, the density peak undergoes 
azimuthal, as well as radial, oscillation about 
the potential minimum. In all but a head-on 
merger (b=0), the displaced gas core does not 
pass directly through the potential minimum, but 
rather undergoes a damped oscillation around the 
potential minimum. Whether due to a head- 
on or off-axis collision, the resulting oscillatory 
motion of the lowest entropy core gas within 
the gravitatio nal poten tial is called “sloshing” 
( Markevitch et all 200ll ). 

Despite the dependence of the long term evo- 
lution of these cold fronts on merger parameters 
(cluster mass ratios, impact parameters, etc.) and 
the plasma properties (e.g. viscosities, magnetic 
fields, etc.), there are several qualities of these 
sloshing cold fronts that appear universal and are 
most relevant to RXJ1347. First is that, based 
on the simulations, there appears to be a mini- 
mum time after subcluster pericenter for the first 
cold fronts to form. This is because the dense 
gas core in the primary cluster must first slingshot 
past the cluster’s potential minimum and will only 
then form the first cold front as it falls back to- 
wards the potential minimum. Though the ability 
to detect the first sloshing cold front formed is lim- 
ited by the simulation spatial resolution, over an 
order of magnitude in both mass ratio and impact 
parameters, simulations show the first cold front 
does not form until «0.3 Gyr after the c losest ap- 
proach of the subcluster (lAscasiba r fc Ma rkevitch 


20061 : ZuHone et al. 2QloT Roedieer et al. " Eoil ) . 
Second, the location of the cold fronts reveal when, 
and from which direction, the sloshing disturbance 
originated. The radial distance (projected or not) 
from the X-ray peak where the cold fronts are 
found always increases with the time since cold 
front formation. That is to say that the sloshing 
cold fronts which formed first are always further 
from the oscillating core than fronts which form 
at a later time. Also, simulations show that the 
first cold front always forms on the opposite side 
of the cluster as the location of the initial distur- 
bance. This is because the initial displacement of 
the primary cluster’s gas core occurs on the same 
side as the pericenter of the disturbing object with 


the first cold front then forming as the core falls 
back from this initial displacement. And finally, 
simulations show that the inclination of the clus- 
ters’ merger axis with respect to our line of sight 
dictates whether the cold fronts will appear to us 
as a spiral inflow (for a merger in the plane of 
the sky), concentric crescent shapes (for a merger 
along the line of sight), or some combination of 
the two resulting from an inclined spiral inflow. 
We discuss the implications of these for the his- 
tory of RXJ1347 in $Q1 


2. Previous Observations of RXJ1347 


As the most X-ray luminous cluster known 
(with Ay =6.0 x 10 45 ergs s ~ 4 from 2-10 keV; 
Voges et al. 1999h . RXJ1347 has been extensively 
studied in the radio, millimeter, s ubmillimeter 
( Gitti et aPl2007a : Mason et al.l20ld ). in the opti- 


cal via both weak and strong lensing (jBradac et al 


2008: Miranda et al 


strong J 
120081) . 


as well as t hrough 


spect r oscopy of memb er galaxies (jCohen fe Kneib 


2002; Lu et al.l 120101) . and of course in X-rays 


( Allen et al. 2002; Gitti et, ahl l2007bl) . Indeed, 


one could argue that of those clusters beyond 
z=0 .3, only the “Bulle t” cluster (IE 0657-56, 


Markevitch et al 


received more at- 
To establish the context and relevance 


(l2004l) f has 

tention. 

of our interpretation, we review here some of the 
previous observations of RXJ1347. We focus pri- 
marily on those contemporary works most relevant 
to our discussion. 


2.1. Previous X-ray Observations of RXJ1347 

RXJ1347 (see Fi g. CD) was first observe d in X- 


rays with ROSAT ([Schindler et alJ 19951 ). With 


the arcminute angular resolution of the ROSAT 
PSPC, RXJ1347 appeared to be a relaxed sys- 
tem with a central cooling peak. Since ROSAT, 
the cluster has been observed by near ly every X- 
ray satellite, in cluding Suzaku dOta et alJ boost) . 
X M M-Newton ( tGjttijdalJ 2007b|h and Chan- 
dra ( Allen et all 2002 ; Bradac et alJ 20081) . Con- 
clusions regarding both the overall gas density, 
temperature and metallicity are in agreement. 
RXJ1347 exhibits a hot (kT>10 keV) ICM, with 
an extension to the sout heast. Prior to the 
Chandra observations ( Allen et al.l 12002 ). there 
was little indication of this southeast structure 
(RXJ1347-SE hereafter) and, excluding this quad- 
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Fig. 1. — Combined 73 ks Chandra image of RXJ 1347. 5- 1145 from OBSIDs 506, 507, and 3592. The image is shown in the 
0.5-2. 5 keV band and has been exposure corrected and background subtracted. The surface brightness edges can be seen on 
either side of the emission peak. The two small circles denote the locations of the primary and secondary cD galaxies (cDl and 
cD2 respectively) and the small blue circle denotes the location of the X-ray peak. In all images, north is up and east is to the 
left. 


rant, the cluster gas distribution appears sym- 
metr ic out to la rge scales (>1 Mpc), leading 
Allen et ah] ( 20021) to conclude that “the cluster 
appears relatively relaxed.” High resolution imag- 
ing and spectral observations, with both XMM- 
Newton and Chandra , reveal a more complex pic- 
ture with several features not typically seen in 
relaxed clusters in hydrostatic equilibrium, includ- 
ing a non-isotropic temperature distribution and 
two extremely hot regions (kT> 20 keV) within 
100 h~ 1 kpc of the X-ray brightness peak (see 

Fig. S). 

2.2. Optical Observations of RXJ1347 

Cohen fc Kneibl ( 2002 ) and Bradac et al.l ( 2008 1 


show that the cluster hosts two cD galaxies, one 
0.5" (3.3 h _1 kpc) NE of the X-ray peak (labeled 
cDl in Fig. []]) and the other, ~ 18" (118 h~ l kpc) 
to the east, which has no detected X-ray counter- 
part (labeled cD2 in Fig.Q]) . The recessional veloc- 


ities for these two galaxies differ by only ~ 200 km 
s -1 ( Cohen fe Kneib 2002 ; Lu et al. 2Q10h which 
is much less than even the lowest estimate for the 
velo city dispersion of the c luster (cq,=910±130 km 


=-i. 


Cohen fc Kneibl [2002 ). 


A dedicated study was performed using 400 
redshifts in the region surrounding RXJ1347 


(ILu et al.1 120101) which, alo n g wit h the previous 
efforts by Cohen fc Kneibl ( 2002f h brought the 
number of cluster me mbers co nfirmed via red- 
shifts to ~140 galaxies. !Lu et al.l (2010) note sev- 
eral new structures in velocity space, located co- 
linearly in projection along a roughly NE /SW line 


through the cluster center (see Fig. 1 in Lu et al 


l20i nh . Among these galaxy overdensities is a sig- 
nificant (7a) peak in the redshift distribution cor- 
responding to a cluster of galaxies approximately 
~7 /i _1 Mpc (projected) SW of the main cluster 
(referred to as RXJ1347-SW hereafter). Although 
this structure appears to be much more distant 
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than RXJ1347 (its mean redshift is ~4000 km s -1 
larger), there is evidence that it is physically as- 
sociated with RXJ1347, either through a previous 
interaction, which is unlikely given the separation, 
or by lying within t he same cosmi c filament as 
RXJ1347 (see §4.2 in lLu et al.llioToi for a discus- 
sion). RXJ1347-SW lies outside both the Chandra 
and XMM-Newto n observations for the cluster, so 
we use ROSAT ( Voges et al.l Il999h to place an 
upper limit on its luminosity of Lx < 7 x 10 44 
ergs s _1 which, using the Ly /M nn ,., r,nn scaling re- 
lations in IZhang et al.1 ( 201 lh . predicts an upper 
limit to its mass (M gaSj 5 oo < 10 14 Mq) t hat is in 
agreem ent with the mass estimate from iLu et all 
( 2010h . Ano ther nearby overden sity of galaxies 
was noted in Bradac et al.l ( 2008h immediately to 
the SW of RXJ1347, coincident with a mass exten- 
sion seen via weak+strong lensing. We also do not 
detect this object in the Chandra image and place 
an upper limit on its luminosity of Lx < 1.3 x 
10 42 ergs s _1 at the distance of RXJ1347. 

There is some disagreement among the different 
lensing analyses as to the presence of the second 
massive dark matter halo that would be associ- 
at ed with eithe r cD2 or t he gaseous subcluster . 
In Bradac et al.l ( 20081) and lMiranda et al~ ( 2008 ). 
their projected mass density maps show an ex- 
tension in the direction of the gaseous subcluster 
30" (^2 70 kpc) southeast of the X-ray bright- 


ness peak. Miranda et al. (2008 ) 


give an approxi- 
mate mass for this structure comparable with the 
primary cluster mass (~1O 15 M 0 ). On the other 
hand, the weak lensing analysis bv lLu et ah ( 2010 ) 
shows no such structure, and does not even show 
an extension of the weak lensing signal in that di- 
rection. We discuss this further in terms of our 
proposed merger history for this cluster (see )j4]). 

2.3. Radio Observations of RXJ1347 


Komatsu et al. ( 19991) and Pointecouteau et al 


(119991 120011) made the first measurements of the 
Sunyaev-Zeldovich effect in RXJ1347. The high 
ang ular resolution ( 20" F WHM at 150 GHz) of 
the [Komatsu et al.1 ( 2001 ) observations showed a 
strong enhancement of the SZ signal 20" south- 
east of the cluster center, at the same location 
as the surface brightness extension observed in X- 
rays (see Fig. [2] a, d). Since the SZ intensity is pro- 
portional to the integrated gas pressure along the 
line of sight, this strong SZ signal was interpreted 


as being due to hot gas, probably shock heated 
by the merger of the subcluster with the primary 
cluster. The presence of hot (kT > 15 keV) gas 
in the so utheast region has been reported from 


Chandra (Allen et al. 2002 


XMM-New ton ( Gittiet^l] 2007b ) and Suzaku ob 


Miranda et al.l 120081) . 


servations ( Ota et al.ll2008 ). and is confirmed here. 

Rec ently jMason et al. ( 2010h and Korngut et al 
( 201ll) reported 90 GHz observations of RXJ1347 
with 10" angular resolution performed with the 
MUSTANG bolometer array on the Green Bank 
Telescope (GBT). In addition to confirming the 
enhanced SZ effect in the region of the subcluster 
merger, the MUSTANG observations (see Fig.[2jl) 
show a pronounced asymmetry in the projected 
pressure, with a high pressure ridge running north 
to south over ~ 0.5' (~197 h~ x kpc) approxi- 
mately 15" (~99 h~ x kpc) east of the X-ray peak. 
The southern end of this ri dge marks th e peak 
of the SZ decrement, which Mason et al.l ( 2010l) 
interpret as a shock between the primary cluster 
and infalling subcluster (RXJ1347-SE). 

We also note that iGitti et ahl ( 2007a ) ob- 
served a radio mini-halo, consistent with the 


sloshing scenario dMazzotta fe Giacintuccil 2008; 


ZuHonc et al] 1201 lh . It is ~420 h 1 kpc in ex- 
tent and elongated to the southeast towards the 
location of RXJ1347-SE, with total flux density 
55.3±0.6 mJy centered on the cluster core. The 
highest resolutio n VL A image (1.7" x 1.2" restor- 
ing beam; Gitti et ah 2007al ) also shows extended 
emission in the core of RXJ1347, with a diffuse 
(i.e. point source subtracted) radio flux density of 
25.2±0.3 mJy. 

3. Our New X-ray Data Analysis of RXJ1347: 
ACIS-I Data Reduction, Imaging and 
Photometry 

RXJ1347 was observed for a total of 80 ks over 
three Chandra observations (Obsids 506, 507, and 
3592). The first two observations were performed 
with ACIS-S and the third with ACIS-I. We re- 
processed the level 1 event files as described in 


Vikhlinin et al 


( 20051) 3 The background flare 
rejection was done using th e ratio of the 2.5-7 
keV and 9.0-12 keV bands ( Hickox fc Markevitch 


1 A11 spectral extraction and spectral analysis were per- 
formed usi n g CIAO v4.3, the CALDB 4.4.1 and Sherpa 
dFreeman et aI1l200ll) . 
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20061) . excluding regions with point sources and 
cluster emission. This resulted in a combined 
~73 ks clean exposure, on which we performed 
the further image processing tasks described be- 
low. Among the products created in this initial 
reduction are background images and exposure 
maps which incorporate the telescope aspect solu- 
tion, vignetting, energy dependent effective area, 
and position and energy dependent detector effi- 
ciency. The background subtracted, exposure cor- 
rected image used for all subsequent position and 
imaging analyses is shown in Fig. [TJ Fig. EKb) 
and Fig. [3] We detected point sources using the 
wavelet detection algo rithm wvdecomp outlined in 


Vikhlinin et al.1 (119981) . The 30 point source re- 
gions within 2.5' of the X-ray peak were checked 
by eye and then excluded from all subsequent anal- 
yses. 


Proper estimation of the X-ray background is 
critical for spectral analysis of the faint extended 
emission from clusters. To measure this, we used 
the blank sky datasets from the Chandra calibra- 
tion database (CALDB v.4.4.1). The aspect solu- 
tion for our observation was applied to the blank 
sky exposures, remapping the background events 
to our source coordinate frame. We then per- 
formed an overall background normalization by 
measuring the ratio of the count rates at high ener- 
gies (9-12 keV, where the telescope’s effective area 
is small) in our source images to that in our back- 
ground images. For each of our spectral extrac- 
tion regions, an identical region on the detector 
was chosen from the blank sky image to provide 
a background estimate. Associated ancillary re- 
sponse files (ARFs) and redistribution matrix files 
(RMFs) were created using CIAO v4.3. 


3.1. Gas Properties Across the Surface 
Brightness Edges 

From the source, background and exposure 
maps, we create surface brightness images with 
different pixel binning: full resolution (no bin- 
ning) and block n with n x n pixel binning (where 
n=2,4,8). The full resolution and block two im- 
ages reveal two sharp edges east and west of the 
X-ray brightness peak (see Figs. [T] and [5Ji). 

In RXJ1347, the presence of two surface bright- 
ness edges to the east and west indicate that the 
core gas may be sloshing, however it is not com- 
pletely clear from which angle we are viewing it. 


We may be viewing an inclined spiral pattern, 
which is turning in an inflowing counterclockwise 
fashion beginning at a radius of ~ 10" (distances 
measured with respect to the X-ray peak) and po- 
sition angle of 140° (measured counter clockwise 
relative to west on the sky) to ~ 5" with posi- 
tion angle 52°. From this perspective, the spiral 
appears somewhat elongated along the N-S direc- 
tion. Such an elongation could indicate either a 
highly anisotropic ICM or gravitational potential 
into which the gas is expanding or, more likely, 
that the perturbing object came along a trajec- 
tory that was inclined with respect to the plane of 
the sky (see 34.31 for a discussion). 

Rather than the spiral pattern, if what we are 
viewing is instead the concentric crescent shaped 
edges associated with a merger more along our line 
of sight, we then face the difficulty in explaining 
the proximity of both cD2 and the gaseous sub- 
cluster to the southeast, which would be very chal- 
lenging indeed. We discuss the observational con- 
straints and implications of both the line of sight 
and plane of the sky scenario in more detail in 

m 

From the surface brightness images, radial sec- 
tors were chosen to encompass the east and west 
edges. We measure the radial surface brightness 
profiles within these sectors and fit them with a 
two component model for the gas density with a 
power law inside the surface brightness edge and 
a /3-model outside with the form: 


(!) 

where rj is the radius of the surface brightness 
edge, 7 i is the power law slope interior to the 
edge, r c is the core radius of the j3 model, ni com- 
bines the model normalization with the emission 
measure normalization for each model in their re- 
spective regions, and J ne is the magnitude of the 
discontinuous density jump at the edge. 

The square of this density model is integrated 
along the line of sight and compared to the 
observed surface brightness profile using a chi- 
squared minimization. Due to the temperature 
dependence on the emissivity of the gas, and thus 
the surface brightness, this is an iterative process, 
where we first fit the surface brightness assuming 
a uniform temperature. Then, upon deprojecting 


p(r) = nl 


for r < r 


T f.njT'i 

J »e V rj ) 

[! + (—) for r > r 
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Fig. 2 . — a) Same image as in Fig. [T] with annular sectors overlaid for surface brightness and spectral extraction b) Same as 
in Fig. [l] with 2x2 pixel binning, showing the spatial extent of the annular regions along with the large scale emission from the 
cluster, c) Temperature map produced as described in £|3.2I d) MUSTANG 90-GHz SZ map (Korngut et al. 2011) showing 
the SZ decrement ridge along the eastern side of RXJ1347. Colors in (d) are of SZ decrement (lower value indicating stronger 
integrated pressure) measured lowest (black) to highest (white) with a logarithmic color stretch. We also highlight the shock 
front which shows as a hot region in the temperature map and an SZ depression (indicating high integrated pressure) in the 
SZ map. All images except (b) are shown at identical positions and scales. Color bar is linear scale in temperature going from 
cool (purple) to hot (white). 


the temperature profile using the resulting density 
model (see H3.3I) . we use the deprojected tempera- 
ture profile to again fit for the emissivity corrected 
surface brightness. Convergence is obtained typi- 
cally after only one such iteration. 

We plot the resulting radial profile for the east- 
ern surface brightness edge (red points, labeled 
“east” in Fig. |4ji) , and find that the surface bright- 
ness profile across this feature may be adequately 
represented by Eq. HI with a jump of J ne =2.2±0.6 
(see Table d at the radius of discontinuity (r ~ 


9.2" ~ 60 h 1 kpc). The density model is plot- 
ted in Fig. 0}:. As discussed in (O given suffi- 
cient time since the initial sloshing disturbance, 
a surface brightness discontinuity on one side of 
the cluster should be accompanied by a second 
discont inuity on the opposite side at a different 
radius ( Ascasibar fc Markevitch 20061 and 31. II in 
this work). Indeed, the western surface bright- 
ness jump lies ~ 4.6" (~ 30 h -1 kpc) from the 
X-ray peak, closer than the eastern edge which 
lies ~ 9.2" (~ 60 h -1 kpc) from the X-ray peak. 
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Fig. 3. — Surface brightness image as in Fig.[T] showing the regions used for estimating the subcluster emission. Each circular 
region is identical in radius, with their centers equidistant from the cluster X-ray peak (shown with vectors pointing from each 
circle). The source region used for the subcluster emission is labeled (blue, solid), whereas the other circles (blue-dashed) are 
used to represent the cluster emission at that same radius. 


Surface brightness and density plots are shown in 
Fig. SJqc (blue points, labeled “west”). That it is 
closer to the X-ray peak indicates that the western 
edge is the younger of the two, resulting from the 
most recent passage of the core gas across the clus- 
ter potential minimum. This feature is also well fit 
by a two component density model, with a power 
law inside the edge, a /3-model outside (Eq.Q]) and 
jump magnitude J Ue = 2.CTto'5- 


3.2. Temperature Maps 

In Figure [2jc we present a temperatu r e ma p 
for RXJ1347, also shown in ZuHone et al. ( 2010h . 
Temperature maps reveal interesting structures 
and inhomogeneities in the cluster gas which 
might otherwise be overlooked by only examin- 
ing the X-ray surface brightness distribution. 


To generate the map , we use the method of 
Markevitch et all ( 20001 ). in which a series of 6 


images in non-overlapping energy bands (0.5-1. 5, 
1.5-2. 3, 2. 3-3. 5, 3.5-4.5, 4.5-6.5, 6.5-9. 5 keV) are 
created along with corresponding exposure maps. 
Each surface brightness image is then smoothed 
identically with a variable length scale Gaussian 
kernel, such that the smoothing scale at each pixel 
decreases roughly as the square root of the number 
of counts in that pixel (out to a brightness of ~0.02 
cts arcsec 2 , below which the smoothing scale is 
constant). This method has the effect of apply- 
ing a smaller smoothing scale in areas of bright 
cluster emission and a larger smoothing scale in 
the low surface brightness regions, resulting in ap- 
proximately the same relative statistical accuracy 
across the interesting bright regions of the cluster. 


Surface Brightness 


Electron Density 



1 10 100 
R [arcsec] 



Pressure 



Fig. 4. — Radial quantities plotted for east and west regions shown in Fig. EJi.b. a) Radial surface brightness, b) deprojected 
temperature, c) electron density, and d) pressure profiles for the regions east (solid x’s, medium (red) shading) and west (open 
dashed squares, dark (blue) shading). Vertical dashed lines are the radii of the edge in each sector. Spectral extraction and 
fitting is described in ED Shading in c) and d) are from the 90% confidence intervals on the model parameters. 


The surface brightness images are then used 
to fit a MEKAL model spectrum for each image 
pixel, binned to those same energy bands, and 
a chi-squared minimization is performed allowing 
only the model temperature and normalization to 
vary. Because this method does not assume any 
a priori shape to the spectral region, as opposed 
to tess ellation or contour binning methods (e.g. 
Sandersl 2006l l. we may be confident that any fea- 


tures shown, particularly near the center, reflect 
an accurate (projected) morphology of the cluster 
gas. 


The temperature map shows large scale inho- 
mogeneities of the cluster gas, indicative of an un- 
relaxed cluster. Of particular note in this map are 
the cluster cool core, the presence of the eastern 
and western cold fronts seen in the surface bright- 
ness image, two very hot (kT>25 keV) regions, 
and the shock front between the primary and sub- 


clusters, also seen in the SZ map (Figure EJi). 


3.3. Spectral Extraction and Fitting 

After fitting the X-ray surface brightness across 
the edges, we extract spectra in these same sectors 
(see Fig. [2] a, b). The spectral bins are necessarily 
wider than the surface brightness bins as we re- 
quire at least 2000 source counts (background sub- 
tracted) in the 0. 8-7.0 keV range per spectral bin 
to constrain a thermal model with small enough 
uncertainties to identify the temperature jumps 
predicted at the cold front edges. Spectral extrac- 
tion and fitting were perfo rmed in an analogou s 
way to our previous work ( Johnson et al. 2Q10h . 
using the CIAO suite of tools. Spectra in each 
spatial region were binned with a minimum of 30 
cts per energy channel and the models were fit 
over the range from 0. 8-7.0 keV. For each spec- 
trum, we fit a thermal MEKAL model convolved 
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with an intergalactic absorption model (xsmekal * 
xszwabs). The cluster’s redshift was fixed at the 
median redshift of its clu ster members (z=0.451) 
determined by Lu et al.1 ( 2010h with an absorbing 
hydrogen column density N h = 4.85 x lO 20 cm -2 
(Dickey & Lockman 1990) at the cluster position. 

We account for the gas lying in projection along 
the line of sight by using a variation of the algo- 
rithms projct for XSpec and deproject for Sherpa. 
Our algorithm takes a series of radially extracted 
spectra, fits a model to the outermost bin first and 
then moves inward. For each inward step, an addi- 
tional model is added to the fit, whose gas temper- 
ature is fixed at the values obtained from the outer 
bins. The difference between deproject , projct and 
our method lies in our separate derivation of the 
3D gas density rather than leaving the density as 
a free parameter in the model fit. The density for 
each radial bin is calculated as described in 
and the ratio of the emission measure in the pro- 
jected region to that in the emission region is used 
as a normalization constraint for each of the pro- 
jected regions. We present these radial quantities 
in Table [1] 

It should be noted that due to the small ef- 
fective area of Chandra at high energies, it is 
difficult to constrain high temperatures (> 20 
keV) with confidence, even with high numbers 
of counts. This is evidenced in the 90% confi- 
dence limits on the temperatures in the hottest 
radial bins (see Table [1]). Despite this limita- 
tion, we do find that the eastern and western 
edges identified in the surface brightness images 
are indeed cold fro nts, as expected from th e ga s 
sloshing model ( Ascasibar fc Markevitchl 120061) . 
From the deprojected temperature and density 
profiles, we are able to calculate the pressure pro- 
file across each edge which, in the absence of high 
radial vel ocities, should be app roximately contin- 
uous ( Ascasibar fe Markevitchl l2006h . In order to 
produce this profile we multiplied the mean tem- 
perature by the mean density in each spectral bin. 
The mean pressure profile does not appear very 
smooth, however this is due primarily to the large 
spectral bins required to constrain the tempera- 
tures, and the large uncertainties in those temper- 
atures. Deeper X-ray observations are necessary 
to decrease the uncertainties in the spatial temper- 
ature distribution. We note, however, that even 
with the coarse temperature resolution, including 


both the temperature and density uncertainties, 
each of the profiles is continuous at the 90% con- 
fidence level (see Fig. 0Ji). 


3.4. The Gaseous Subcluster 


The gaseous subcluster to the southeast of the 
X-ray core lies within the halo of RXJ1347 (see 
Fig. HD, making a study of its spectral proper- 
ties problematic. We isolate its emission by se- 
lecting a circular region in the surface brightness 
image containing its emission (see Fig. H along 
with three other identically sized regions at the 
same distance from, but at different position an- 
gles with respect to, the X-ray peak. The flux 
we measure in the region surrounding the sub- 
cluster includes some cluster emission, which we 
measure as f x , su b+ciust = 3.8 ± 0.5 x 10 -13 ergs 
s^ 1 cm -2 and compare with the mean flux in the 
other three regions (j x , other = 2.3 ± 0.4 x 10~ 13 
ergs s _1 cm -2 jl. So our estimate of the sub- 
cluster flux is the difference of these, or f XtSU b = 
1.5 ± 0.9 x 10” 13 ergs s _1 cm -2 . This subcluster 
flux translates to a subcluster luminosity (at the 
distance of RXJ1347) of L XjSub = 1.1 ± 0.6 x 10 44 
ergs s" 1 . Us ing the L-M gas scaling relations from 
Zhang et ah d20Tl|j 2 3 l . this corresponds to a sub- 


cluster gas mass of M goSjSu{ , ~ 4.8 ± 2.5 x 1O 13 M 0 
within rs nn- Assuming a core gas mass fraction 


f gas = 0.1 ( Vikhlinin et al. 2006;), this gives a mass 


for the subcluster of M 50 o,sub = 4.8±2.5x 10 14 M Q . 
Given the merger scenario we propose in (J4l the 
subcluster would have likely lost some portion of 
its mass during its first pericentric core passage. 
This suggests that our mass estimate is probably 
a lower limit. 


4. Discussion - The Merger History of 
RXJ1347 


While previous works agree that RXJ1347 ap- 
pears to have undergone a sign ificant merger 


( Gitti et al. 2007a : Lu et al. 2010h . there is lit- 


tle consensus on the details. The cluster has a 
complex morphology at nearly every wavelength, 
and there is no single merger scenario thus far 


2 A11 fluxes here computed in the 0.5-2. 5 keV energy band. 
Errors are systematic variation of la based on the three 
regions surrounding the cluster 

3 There are separate scaling relations for both disturbed and 
undisturbed systems, and we use the former. 
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Table 1 


Radial Temperatures for RXJ1347 


Region 

R 

kT 

X 2 (dof) 


["] 

[keV] 


West 

70.0 

15.4+2:7 

191.7 (197) 


25.8 

24.6 ±11 

230.4 (297) 


9.5 

24.0±|'f 

195.2 (257) 


2.5 


87.0 (108) 

East 

95.2 

i5.o±r 4 

94.0 (80) 


34.9 

26.8+g 1 j 6 

110.7 (136) 


13.7 

12.5±fy 

191.2 (214) 


4.5 

0 1 +1-3 

O’-L-l.O 

112.9 (160) 


Note. — MEKAL model temperatures for 
the regions shown in Fig. [2] (a,b). (col 1) 
Spectral extraction region identifier, (col 2) 
Emission weighted radius of spectral bin. 
(col 3) Deprojected MEKAL model tempera- 
ture. Errors on kT are the 90% confidence in- 
terval, A% 2 =2.71, for one interesting param- 
eter). (col 4) X 2 f° r model and number of 
degrees of freedom. 


Table 2 


Edge Properties for RXJ1347 


Cluster ID 

Edge Radius ["] 

T 2 

Tl 

^el 

n e2 

Pi 

p 2 

East 

a 9+°- 5 

y '^- 0.5 

1 eo + 0.53 
1,0 - 0.40 

c\ qq -|- 0.60 

z.zz-o 58 

1 40 -1-0 ' 65 
i ‘^ U - 0.50 

West 

4 5+ 02 

3 24 +2 ' 41 

1 Q 8+ 0,61 
1 - yo - 0.45 

0 6 i + 0 - 54 
uo ± — 0 .28 


Note. — Values in each column are the ratios of various gas quantities 
(Temperature, gas density and gas pressure) across the eastern and west- 
ern surface brightness edges in RXJ1347 (see Fig. [Til, where Ni represents 
the value of the gas property N inside the front (r < r ec ig e ) and N 2 is the 
value outside the front (r > r e d ge )• The errors presented on the ratios 
are based on the 90% confidence limits from those quantities. 
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presented which can explain the unique charac- 
teristics of this cluster. Specifically, the key ob- 
servables that must be explained by an acceptable 
merger model are: 

• The presence of two sloshing edges: one at 
~ 9.2" (^60 h~ l kpc) east of the X-ray peak, 
and a second ~ 4.6" (~30 h~ x kpc) to the 
west of the X-ray peak. 

• A very hot region located ~ 17" (~113 h~ l 
kpc) to the SE of the primary cluster. 

• A second cD galaxy with no detected X- 
ray gas atmosphere located ~ 18" (118 h” 1 
kpc) directly east of the primary cluster and 
~ 11" (72 h~ l kpc) north of the gaseous 
subcluster. 


4.1. Previous Merger Scenario 


The presence of the high pressure ridge along 
the east side of the cluster (see Fig.[2ji), connecting 
the second BCG (cD2 in Fig. |TJ) to the gaseous sub- 
cluster, sug gests that these obj ects are all related. 
As noted in Mason et al.1 (120101) , these features are 


both consistent with cD2 having moved north (in 
projection) along the east side of RXJ1347, with 
the gaseous subcluster having been either stripped 
or shock heated (or both) as a result of this pas- 
sage. This is the simplest model that agreed, un- 
til now, with all of the observations. One problem 
with this model is that sloshing edges are produced 
well after the passage of the subcluster, which does 
not agree with the current location of either cD2 
or the gaseous subcluster. 


4.2. Sloshing Constraints on the Merger 
History - A Second Crossing 

As discussed in an simulations of gas slosh- 
ing predict a minimum time («0.3 Gyr) should 
elapse between the first pericentric crossing (t per ) 
of the perturbing cluster and the formation of the 
first sloshing cold front, irrespective of the merger 
axis inclination with respect to the line of sight. 
For the subclusteJl to have perturbed RXJ1347, 
it should now be a minimum distance away from 


4 Here we make special distinction between the gaseous sub- 
cluster, RXJ1347-SE, and the subcluster in general, com- 
prising cD2, the gaseous subcluster, and their combined 
dark matter halo. 


pericenter for us to see sloshing cold fronts in the 
core. For example, if we assume a reasonable 
value for the subcluster’s average velocity since 
tper, say 1000 km s _1 , then it should be at least 
~300 kpc away in 0.3 Gyr. This minimum dis- 
tance is further constrained by the observation of 
a second sloshing cold front, which occurs only 
after the density peak has fallen past the poten- 
tial minimum, reached apocenter, and has then 
fallen back againjf] Simulations show that this 
second front forms, at the earliest, ~0.6 Gyr af- 
ter t per, increasing the predicted separation dis- 
tance between pericenter and the current location 
of the subcluster to ~600 kpc. This is consider- 
ably larger than the observed (projected) separa- 
tion distance (~161 kpc) between RXJ1347 and 
the location of cD2 (which should also be the lo- 
cation of the subcluster, as both the dark matter 
and the cD galaxy pass essentially collisionlessly 
through the envelope of RXJ1347). This discrep- 
ancy indicates that the simple scenario where the 
subcluster is merging for the first time is not suffi- 
cient. We therefore propose a new scenario, which 
we will refer to as the “second crossing” scenario. 

4.3. A New Merger History For RXJ1347 

The simple interpretation of RXJ1347 now ex- 
periencing the first crossing of the subcluster does 
not correctly predict the presence of the slosh- 
ing cold fronts observed in the cluster core. It 
is also inconsistent with the simulation result that 
the primary cluster’s core gas does not separate 
from its potential minimum until after the first 
passage of the subcluster. We propose instead 
that the subcluster first passed on the west side 
of RXJ1347, moving in a (projected) southerly di- 
rection, initiated core gas sloshing upon reaching 
pericenter in its orbit, and has now returned for a 
second pass, this time from the south, and much 
closer to the core. During this second passage, 
the subcluster had its core gas stripped as it ap- 
proached pericenter, and is now near this point 
in its orbit. This “second crossing” scenario ad- 
dresses each of the three cluster observables noted 
in (21 

The presence and location of two sloshing cold 
fronts around the primary cluster’s core: By al- 
lowing the core gas sloshing to be initiated on 

5 This is assuming a near plane of the sky viewing angle. 
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the first subcluster passage, we solve the discrep- 
ancy between the minimum time since percen- 
ter (~0.6 Gyr) and the current location of the 
subcluster. In this scenario, we require that the 
time between pericentric crossings of the subclus- 
ter be larger than 0.6 Gyr, as sim ulation s pre- 
di ct i s th e case (c.f . Tittlev & Henriksenl 2005 


i itticy 


Ascasibar fc Markevitc 


ZuHone et all 2010; 


Roediger et al .1120 111 and Hl.ll in this work). Simu- 
lations also predict that the first cold front should 
be produced on the opposite side of the cluster 
from the perturber’s pericentric location. For 
RXJ1347, we propose that the initial subcluster 
pericentric crossing was on the west side of the 
cluster. Thus the first cold front should be on 
the east side. Our measurements of the cold front 
radii (see Table [2j indicate that the eastern cold 
front is farther from the X-ray peak (9.2") than is 
the western cold front (4.6"), e.g. it was formed 
earlier, in agreement with the “second crossing” 
scenario. We also note that the lowest mass sub- 
cluster that has been si mulated to produce s losh- 
ing is of order 10 13 M 0 ( Roediger et all [201 il l. As 
shown in H3.41 our estimate for the subcluster gas 
mass alone is M gas ~ 10 13 Mg, putting its total 
mass Mt ~ 10 14 M 0 . This implies the primary 
cluster is of order ten times the mass of the sub- 
cluster, which is well above the mass ratio limi t 
required to initiate sloshing ( ZuHone et al. 2010h . 

A very hot region located between the primary 
and sub- cluster, and the current location of cD2 
and the gaseous subcluster: Shocked gas between 
the primary cluster and the gaseous subcluster 
is predicted in all merger simulations for which 
the subcluster possesses gas of its own. This 
is because, being initially bound to the subclus- 
ter, the subcluster gas velocity exceeds the sound 
speed of the main cluster as it approaches peri- 
center in its orbit. The subcluster undergoing a 
second crossing will not change this. The shock 
from the subcluster’s first crossing would have 
passed through the primary cluster’s core already, 
likely contributing to the initial displacement of 
the R XJ1347’s gas core fro m the potential mini- 
mum (jChurazov et al .1120031) . Before the subclus- 
ter reached apocenter on its first pass, this shock 
would have traveled well beyond the cluster core, 
so that there would be little evidence of the front 
during the subcluster’s second crossing. 

As the subcluster returned for its second cross- 


ing, this time from the south, its trajectory 
brought it closer to the core of RXJ1347, where 
it is now near pericenter, at ~ 18" (118 h~ x kpc) 
east of the primary cluster center. At this closer 
pericentric distance, the combination of the sub- 
cluster’s high velocity and the ICM density of 
RXJ1347 increased the ram pressure on the sub- 
cluster’s core gas enough to strip the gas from the 
dark matter halo. The gaseous subcluster is still 
traveling supersonically so we see a s hock f ront be- 
tween it and the core of RXJ1347. iBradac et al. 
( 20081 see their Fig. 4) produce the highest reso- 
lution lensing map of the region around RXJ1347 
and show that the projected mass density is elon- 
gated towards the region in between cD2 and the 
gaseous subcluster. That the mass density elonga- 
tion does not lie directly on the gaseous subcluster, 
or on cD2, indicates that some fraction of the pro- 
jected mass density is lying north of the subcluster 
in the direction of cD2. This supports the notion 
cD2 has traveled along with the dark matter halo 
to the north of the gaseous subcluster. 


4.4. Simulated Merger 


A merger scenario similar to the one described 
above was modeled using the FLASH hydrody- 
namic code ( Frvxell et alJl200(l i. The model ini- 
tial conditions are as described in ZuHonel ( 201lh 
for a different set of simulations with mass ratio 
R=10 and impact parameter b~l Mpc. When 
viewing the simulation in the XY plane (the plane 
of the clusters’ merger axis), the subcluster ap- 
proaches on the western side of the main cluster 
on a southerly trajectory. Near pericenter, it ini- 
tiates sloshing in the primary cluster core. The 
subcluster reaches apocenter, then falls back to- 
wards the main cluster. Its proximity to the pri- 
mary cluster’s core coupled with its relative veloc- 
ity causes ram pressure on its leading edge that is 
strong enough to remove the subcluster gas core 
from the subcluster potential minimum. 


In Fig. [5] we show one particular snapshot in 
time of the simulated projected surface bright- 
ness, temperature and mass density of the cluster 
and subcluster and compare these with observa- 
tions of RXJ1347. Given the variety of distinct 
characteristics in RXJ1347 with which to com- 
pare, the qualitative similarities between the sim- 
ulations and observations are remarkable. In both 
the simulated and observed surface brightness irn- 
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100 kpclM.n 




Fig. 5. — Figure comparing a simulated cluster merger to observations of RXJ1347. On the left are the simulation snapshots 
with (a) projected surface brightness, (b) spectroscopic-like projected temperature and (c) projected mass density with linearly 
spaced contours of constant mass density. The observations on the right are (d) 0. 5-2.0 keV surface brightness image smooth 
with Gaussian (FWHM=0.5 // ), (e) temperature map as in Fig. [2] (f) full resolution Chandra surface brightness image showing 
the location of cD galaxies. The spatial scales for each simulated image are the same, as are the spatial scales for each Chandra 
image. 


ages (Fig. [5] a,d) the sloshing edges are apparent 
to the east and west of the primary cluster. The 
distances of the simulated fronts from the X-ray 
peak do appear different than those in the X-ray 
image. Specifically, in the simulated image, the 


eastern edge is closer to the X-ray peak, indicat- 
ing that the sloshing there is at a more advanced 
stage than in the observation. This is not surpris- 
ing given that the current location of the density 
peak is dependent on many factors (pericentric 
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crossing time, impact parameter, mass ratio, the 
subcluster’s relative velocity and orbital trajectory 
to name a few). The subcluster core also appears 
more concentrated in the simulation image com- 
pared to the Chandra image. Other merger re- 
lated structures are visible when comparing the 
projected temperature maps (Fig. 0b, e). Despite 
the poorer resolution of our temperature map rel- 
ative to the simulation, we can still identify the 
location of the sloshing edges. We are also able to 
identify the cool region of the infalling subcluster 
as well as the characteristic shock front heating 
the gas in between the primary cluster and the 
subclu ster. In our temperature map (and i n other 


works; Ota et all l2008t iMason et al. 201Clh . there 


are two very hot (> 20 keV) regions which we do 
not see in the simulated temperature map. These 
hot regions may have been shock heated during 
the first passage of the subcluster. In the simula- 
tions, the entire cluster is heated as a result of this 
initial shock front, and by the second crossing the 
overall cluster temperature has been raised from 
its pre-shocked state. The higher central temper- 
ature in the simulated cluster core is evidence of 
this heating. The two hot regions, not associated 
with SZ decrements, which we find in our temper- 
ature map may indicate that the gas in these re- 
gions has cooled less relative to the simulation (i.e. 
less mixing), or that there has been some other 
mechanism (such as another shock front) which 
has maintained it at such a high temperature. 

Finally, we compare the simulation projected 
mass density to the locations of the two cD galax- 
ies (Fig. 0c,f). Because the galaxies and the dark 
matter in a merger should pass through nearly 
collisionlessly, the distribution of galaxies should 
trace the distribution of dark matter. In particu- 
lar, since a cD galaxy resides at the bottom of its 
cluster’s potential well, it should trace the location 
of the dark matter density peak. Comparing the 
images, we see that the location of cD2 relative to 
cDl coincides with the location of the subcluster 
density peak relative to the main cluster density 
peak in the simulation. The location of the sub- 
cluster halo in the simulation also agrees w ith the 


elongation found in lens ing maps ((Miranda et al 
2008; Bradac et al.ll2008h . 


5. Summary 

We have presented a new scenario for the 
merger history in the most luminous X-ray clus- 
ter, RXJ1347. This scenario takes advantage of 
the wealth of multiwavelength observations for the 
cluster. We perform a new analysis on the exist- 
ing X-ray data that shows the core gas is sloshing, 
with sloshing cold fronts present to both the east 
and west of the cluster core. The morphology of 
the sloshing core gas places additional constraints 
on the merger history of the cluster: namely that 
the interaction which initiated the gas sloshing 
occurred at least 0.6 Gyr ago, and that the per- 
turbing object first passed along the western side 
of the cluster. In order to reconcile this with the 
current positions of the nearby gaseous subcluster 
and second cD galaxy, the perturbing subcluster 
should have first passed along the western side of 
the primary cluster, in a southerly direction and 
initiated sloshing near pericenter in its orbit. We 
are now seeing the subcluster when it is again 
near pericenter, but this time it has crossed close 
enough to the primary cluster’s core that its gas 
has been ram pressure stripped. We make qual- 
itative comparisons between the multiwavelength 
observations for the cluster and high resolution 
hydrodynamic simulations. We find that many of 
the unique merger structures around the cluster 
are well reproduced in the simulations. This par- 
ticularly well studied cluster serves as an example 
of how core gas sloshing may be used to constrain 
the merger histories of galaxy clusters. 
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